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Abstract 



For quantum field theories that flow between ultraviolet and infrared fixed points, central 
functions, defined from two-point correlators of the stress tensor and conserved currents, in- 
terpolate between central charges of the UV and IR critical theories. We develop techniques 
that allow one to calculate the flows of the central charges and that of the Euler trace anomaly 
coefficient in a general N=l supersymmetric gauge theory. Exact, explicit formulas for SU(N C ) 
gauge theories in the conformal window are given and analysed. The Euler anomaly coefficient 
always satisfies the inequality ayy — am > 0. This is new evidence in strongly coupled theories 
that this quantity satisfies a four-dimensional analogue of the c-theorem, supporting the idea 
of irreversibility of the RG flow. Various other implications are discussed. 
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1 Introduction 



In two dimensions there are many known examples of quantum field theories that flow under 
the renormalization group (RG) from a conformal fixed point in the ultraviolet to another fixed 
point in the infrared. Often, it is possible to work out exact results, using properties and 
techniques that are special to two dimensions. The operator product expansions of two stress 
tensors or two conserved currents contain central charges c and k which encode fundamental 
properties of the conformal theories that appear at the limits of the flow. The Zamolodchikov 
c-theorem [|J states the important inequalities cjjv — cm, > and kuv — km > that place 
constraints on the flow and have a useful physical interpretation in terms of RG irreversibility 
and the thinning out of degrees of freedom as one moves to longer distances. 

In four dimensions it is more difficult to establish the existence of conformal fixed points. 
When they exist, a quantum field theory can be described as a radiative interpolation between 
pairs of four-dimensional conformal field theories. The problem is then to identify relevant 
physical quantities and study their renormalization group flow from one fixed point to the other. 
We call this problem the RG interpolation. While the RG interpolation seems to be very difficult 
in the general case, it simplifies considerably for supersymmetric theories, where many examples 
of interacting conformal fixed points have been studied. In particular, supersymmetric gauge 
theories in the "conformal window" ^ have nontrivial IR fixed points; this is also the subset 
of theories with electric-magnetic duality. The general relation between the trace anomaly and 
the chiral anomaly of the .R-current in supersymmetry allows us to solve the RG interpolation 
problem for supersymmetric gauge theories, and the principal application of these techniques 
is to theories in the conformal window. 

Another difficulty of four dimensions is that no analogue of the c-theorem has been proved, 
although the coefficient a of the Euler density in the curved space trace anomaly has been 
proposed § as a c-theorem candidate, and the flow ajjv — a iR is known to be positive in all 
cases where it can be tested. We shall have more to say about this later. 

In two-dimensional conformal theories the operator product expansions of two stress tensors 
T^ v or conserved currents J M are closed, namely no new operators appear, but the situation is 
more complicated in higher dimensions. As in two dimensions, the TT and J J OPE's define 
primary central charges c and b, respectively. But these OPE's are not closed ||, ||. New 
operators X with anomalous dimension appear, and the OPE's of the new operators define 
secondary central charges d and b' ||. AS operator can be used to deform the theory off- 
criticality and its anomalous dimension h then coincides with the critical value of the slope of 
the /3-function ||. This property indicates that information about the off-critical theory can 
be obtained by studying its critical limit. 

It is reasonable to suggest that the central charges c, c', b, b', the Euler coefficient a and 
the anomalous dimension h are the fundamental parameters of four-dimensional conformal 
theories and to study these quantities, which depend on the dimensionless couplings of the 
theory. The lowest order, two-loop radiative corrections to the central charge c were obtained 
by several authors and in particular by Jack || in the most general renormalizable theory of 
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scalar, spinor, and gauge fields. In |5J the secondary central charge c' was computed to two-loop 
order for general N = 1 SUSY gauge theories and Jack's result for c was specialized to the case 
of supersymmetric couplings. The general N = 1 theory contains the extended N = 4 theory 
and the N = 2 theory with the critical number of hypermultiplets as special cases, and it was 
observed that c and c' are constant on the well-known marginal fixed lines of those theories. 
It is a consequence of the c-theorem in two dimensions that the central charge c is constant 
along lines of marginal deformation. This suggests that the central charges are invariants of 
super confer mal field theory in four dimensions (SCFT4) - they do not change within families of 
continuously connected theories. The anomalous dimension h does vary along marginal lines. 

In ref. (8| non-perturbative definitions of primary and secondary central functions were 
given in terms of two- and four-point correlators of the stress tensor and conserved currents. 
The central functions interpolate between the critical values of the central charges. One goal 
of the RG interpolation is to work out nonperturbative expressions for these functions. We 
are going to show that this can be achieved in supersymmetric theories for all primary central 
functions, as well as a special subclass of secondary central functions. 

In particular, the major results of our analysis are non-perturbative formulas for the flow of 
the primary central functions b, c and a, that agree with perturbative calculations and can be 
applied within the conformal window of SU(N C ) theories with Nf flavors to give exact formulas 
for the total flows bjjv ~ bin, cjjy — cm and ajjy ~ a lR i n terms of N c and Nf and no other 
parameters. The first of these is strictly negative so that no 6-theorem holds, while the second 
changes sign from positive to negative as Nf increases from 3N c /2 to 3N C . On the other hand, 
the flow of a is positive in the entire conformal window. A byproduct of our general analysis is 
the marginal constancy of the primary central functions to all orders in perturbation theory. 

At the perturbative level, while the two-loop contribution to the primary flavor central 
charge b turns out to vanish along marginal lines, in agreement with the nonperturbative for- 
mulas, the radiative corrections to the secondary flavor central charge b' are not marginally 
constant. At the moment, we do not have a general nonperturbative treatment of secondary 
central functions. We hope to discuss our two-loop calculations, which are based on an interest- 
ing application of conformal symmetry to calculations of Feynman diagrams, and the situation 
of secondary central charges, in a subsequent paper ||. 

The flow of the flavor central charge b is analysed by two methods, one in Section 2 and the 
other in Section 3. The second method is more direct, but only the first can be easily extended 
to the gravitational central charges c and a as described in Section 4. In Section 5, our formulae 
for the total flows bjjv — bin, cjjv — cir and ajjv — &IR are discussed, and we give there our 
conclusions and outlook. In Appendix A, the component gravitational anomalies are obtained 
from their curved superspace counterparts. 
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2 Flow of the flavor central charge 



We consider SU(N C ) supersymmetric QCD with Nf quark flavors. The theory contains gauge 
superfields V a (x, 9, 9), a = 1, ... , JVj? — 1, whose physical components are the gauge potentials 
A^(x) and Majorana gauginos A a (x). There are also chiral quark and anti-quark superfields, 
Q ai (x,9) and Q a i(x,9), respectively, where a = 1,...,N C and i = l,...,Nf. The matter 
components of Q m are the complex scalars <f> al and Majorana spinors ij) a% , while Q a , L contains 
4> a i and ifiai. This electric theory has the usual gauge interactions and no superpotential. Later 
we will extend the treatment to the magnetic theory where there are additional gauge neutral 
chiral superfields and a cubic superpotential. 

The theory has an anomaly-free SU(Nf)n x SU(Nr)~ x U(1)b global symmetry group, 
and there are conserved currents that appear as the 99 components of the superfields Qt A Q, 
Qt A Q, and B = (QQ - QQ)/N C . Here t A and t A are matrix generators of the fundamental and 
anti-fundamental representations of SU(Nf), respectively. The component Noether currents of 
the Euclidean signature theory are (see for details of the notation) 

J^-^Lt^ + ^D^t^, (2.1) 



Jfl ~N r . 



1 _ 1 ~ - ~ <-> ~ 

2^75^ - -^Ipl^ + <t> Dfj, cj) - <j> (p 



There are also classically conserved but anomalous Konishi and R currents that are given by 

(2.2) 



2^7^75^ + 2^7/^75^ + <t> D fl <p + (f) 



1- 1- ~ ~ 2 - ^ 

= 2 A "7M75A a - -(V7 M 75^ + V-7m75^) + -(<f> <f> + <f> 0). 

The first is the 99 component of K = QQ + QQ, while the second is the lowest component of the 
supercurrent superfield J a a that also contains the stress tensor and supersymmetry currents. It 
is well known that and have internal anomalies Q, expressed by the operator equations 
(see 0) 

The anomaly-free, RG invariant, combination AHI of and R^, namely 




S„ = R„ + - 1-^-7 \K ll , (2.4) 



1 We distinguish between internal anomalies, involving the quantum gauge field, and external anomalies, 
involving an external classical gauge field. 
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will be important for us. The coefficient of is the numerator of the exact NSVZ |TT| (3 
-function 

lb7r z 1 — g z Jy c /6ir A 

and 7/2 is the anomalous dimension of the superfield Q (or Q). 

We now let J^x) denote any one of the conserved flavor currents of ( |2.1| ). Its two-point 
current correlation function has the form 

(Mx)J v (Q)) = ^{d,d u - □ V) ( b{9i ll x)) ) , (2-6) 

Since J^(x) is conserved and thus has no anomalous dimension, the Callan- Symanzik equation 
requires that b(g(l/x)) depend only on the running coupling g(l/x) that satisfies 

^-0W». (2.7) 

As discussed in || b(g(l/x)) is a primary central function that interpolates between flavor 
central charges at the fixed points guv and gm of the renormalization group flow. Specifically 

b uv = lim b(g{l/x)) = b{g uv ) (2.8) 

x — >0 

b IR = Jim b(g(l/x)) = b(g IR ) 



X 



The principal result of this section will be an exact non-perturbative formula for buv ~ bin- 
This will take some discussion, but the basic ideas are simple and we will list them here before 
we begin the derivation: 

1. In the presence of an external source for the current J M , the trace anomaly of the 
theory includes the familiar internal anomaly plus an external anomaly of similar form, 

where F" and B^ u = d^B u — d u B^ are internal and external field strengths, respectively. 
The general form of ( |2.9| ) follows uniquely from gauge invariance, locality dimensional 
analysis, and parity Jl3[ . As discussed in refs. |12|, |l^], the coefficient of (F^ u ) 2 is the 
numerator of the /3-function (|2.5|) . The coefficient of (B^) 2 is the subject of our investi- 
gation. Following the discussion of H, that we review below, q can be identified with 
the coefficient of a local term that appears in the explicit scale derivative of the regulated 
correlator (^lj|). This coefficient is a scale dependent function b(g(^)) that can be shown 
to have the same UV and IR limits as b{g{^)). 



2 As discussed in Section 3, q can be actually identified with the limit as e — > of the U(l) /3-function f3 e (g,e) 
of an extended SU(Nf) x U(l) gauge theory with U(l) coupling constant e, and this provides an alternative 
method to derive the flow. We do not use this observation here, because only the present method can be easily 
extended to the gravitational central charges as discussed in Section 4. 
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The introduction of external sources can be done supersymmetrically by embedding BAx) 
in an external gauge superfield B(x,6,6) coupled to the superspace flavor current. The 
usual relation between the trace and d^R^ anomalies implies that 



In this way the flow of the central function b(g(fi)) is related to the flow of the anomalous 
correlator (R^(x)J u (y)J p (z)) = sS^)dW^)^ R ^ x ))- 

3. The final ingredient is 't Hooft anomaly matching |1J] for the internal anomaly- free current 
in fl2,4| ). Since S p is quantum-conserved in the absence of sources, its external anomaly 
is scale-independent. This implies that the particular combination of the external dR and 



dK anomalies in (2.4) is independent of scale, and a useful non-perturbative expression 



for the flow of b(g{[i)) emerges from this observation. 

To begin discussion of the first, and probably least familiar, ingredient, we consider the 
structure of the correlator (|2.6|) to any finite order of perturbation theory and construct a line 



of argument from which we extract an all-order result. We consider the computation of (2.6) 
in two stages. In the first stage we work at separated points and regulate all sub-divergences 
at the scale \i. To any finite order, b(g(l/x)) can be expressed in the form 

b(g(l/x)) = bn(9(v))t n , t = In(xM), (2.11) 

n>0 



where b n (g) is a polynomial in g. Evaluating ( p. 11 ) at x = l//x, we see that b(g) = bo(g). The 
Callan-Symanzik equations imply 

P(g)b' n (g) + (n + l)b n+1 (g) = (2.12) 

so that all b n (g) for n > 1 are proportional to /3(g) and can be expressed in terms of bo(g), /?(<?) 
and derivatives. 

When Q2.ll ) is inserted in ( |2.6j ) one finds t n /x 4 terms which are too singular to have a 



finite Fourier transform. To correct this we regulate this overall divergence at x = using the 
generalized differential identity || 

(Ins//)" n! " 2^+1 

= "2^1 ° £ + " anKxl (2 - 13) 

which is an application of the method of differential renormalization fl5|| . Any other method 
in which the overall divergence can be separated from sub-divergences would also work. For 
example, dimensional regularization in x-space might be used. The right hand side of ( |2.13 ) 



contains the Lorentz invariant solution of the differential equation Df(x) = t n /x 4 which is 
unique up to the additive numerical constant a n . Combining ( 2.1 1] ) , Q2.13|) and (pTq), we find 
the fully regulated correlator 



n ! " 2 k t k+1 1 



(2.14) 
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The explicit scale derivative of ( ^6| ) gives |T(| the correlator (J M (x)J y (0) / d 4 z@(z)) where 
O is the trace anomaly Q2.S| ) . We compute this scale derivative as \i d/dfi acting on |U4|). The 
result is expressed as the sum of the local contribution of the k = term plus a non-local term 
proportional to /3(g) because of ( |2-12j ) 



d b(9{lM)_ = 27V 2~ b{gifjL))§{x) + p {g0l))n *M ^ (2 . 15) 

CJ LL X X 



where 

— m _ 77 I 

b(g(n)) = Y,K(9(»))^ ■ (2.16) 

n Z 

F(x) contains the sum of all k > 1 terms in the scale derivative, and the a n have dropped out 
because they have vanishing scale derivative. The decomposition between local and nonlocal 
contributions in the above expression is not universal. For example, one can add an arbitrary 
x -independent function A(g) to F(x) and redefine b(g) as b(g) — (3(g)A(g) in the local part. 
However, such contributions are proportional to (3(g(fj,)) and vanish at the fixed points, so the 
total flow of the central function b(g(l/x)) can be computed from the flow of the function b as 
defined in Q2.16Q . 

It is easy to derive, using ( 2,12 ), a differential equation for b(g), namely 

P(g)^- + 2b(g) = 2b(g), (2.17) 



which is solved by 

2 

%(M)) = ~ 2 f 1 dfyg^)). (2.18) 
V z Jo 

These formulas have been derived within perturbation theory, but we shall regard them as non- 
perturbative results. Either formula shows that the functions b(g(/j,)) and b(g(/j,)) coincide at 
fixed points of the RG flow, as well as to the second loop order in perturbation theory around 
the free fixed point 

Finally, one can identify b(g(/j,)) with the coefficient q appearing in eq. ([2.9|). To show this 
we write the generating functional for the current correlation functions as the schematic path 
integral 

e -r[B M ] = j [d§] e S[*]-H J d+xJ^B^x)^ ( 219 ) 

The source couples just like an abelian gauge field without kinetic term. The scale derivative 
jjL-g- corresponds to the insertion of / d i zQ(z) inside the path integral, so that 



M |_ e -r = J me -sm + ifd*xMx)B» { x) J d 4 z 



SN c -N f (l-j) 1 

^2 y p ^) + ^( B ^) 



(2.20) 



To two-loop order there are no lnx[i terms in (2.11 ) so that b(g) = b(g). 



7 



with similar internal and external contributions. From ( |2.20| ) we see that the scale derivative 
of the current correlator satisfies 

d_ 



^— ( J p (x)J u (0)) = (J p (x)J u (0) I d A zQ(z)) (2.21) 



q{d,d v - U8^)8\x) - 3Nc ^ 7) (J,(x)M0) J d*z{F° v f). 



It has the same form as ( |2.15 ), namely the sum of a local term plus a non-local term proportional 



to (3(g(/j,)), so that we can identify the coefficients of the local terms, again up to contributions 
0(/?(g(/u)) that vanish at fixed points. (Note that the last correlator in ( 2.21| ) is 0(g A (/i)) and 
vanishes in the ultraviolet.) We can thus write 

q=^b(9(»)), (2-22) 

with the understanding that possible 0((3(g(fi))) corrections, that are irrelevant for the total 
flow of 6(<?(/i)), are omitted. 

The next step is to use the d p R^ anomaly ( 2.1 0| ) which is the supersymmetric partner of 



the trace anomaly (|2.9[ ) to compute the flow of b(g{[i)). We need the anomalous correlation 
functions (S p (x)J u (y)J p (z)), (R p (x) J u (y) J p (z)) and (K p (x) J u (y) J p (z)) and their flows. The 
anomalous divergence of such current correlators is one-loop exact only if all currents have no 
internal anomalies. Otherwise the so-called "rescattering graphs," which contain an internal 
triangle in which one of the currents communicates to a pair of gluons, are responsible for higher 
order non-local corrections [Oj. Therefore 't Hooft anomaly matching holds for the correlator 



(SJJ) where S p is the anomaly- free current of (2^), but not for (RJJ), and in general not 
for (KJJ) . The anomalous Ward identities of these correlators can be written as the following 
equations for matrix elements of S p , R p and K p in the presence of the current source B p 

(d p S p ) = ^s B pu B p „, (2.23) 
(d^R/j,) = -^^b(g(n))B pi/ B pi/ + ■■■, 

(d p K p ) = --L^k(g(^))B pu B pu + ■■■. 

Here sq is independent of scale, while the d^K^ anomaly coefficient is defined as the scale- 
dependent function k(g(fi)). In the (d p R p ) equation "+..." indicates the non-local contribution 
of the internal anomaly term of (|2.10 ), There is a similar non-local contribution to (d^K^), 



which cancels that of (d p R p ) in the linear combination (^J) that gives (d^S^). These contribu- 
tions are irrelevant for our analysis, since the possible local terms they contain are 0(fl(g(n))). 
From (|2.4| ) we see that the local terms in ( 2.23| ) satisfy 



b(9)+ (1-^-7(5))%) = so- (2.24) 
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In applications to asymptotically-free electric supersymmetric QCD, the coupling g(n) vanishes 
in the ultraviolet, and so does 7(5); thus the ultraviolet contributions buv and kuv can be easily 
obtained from 1-loop contributions to the relevant 3-point correlators. These are normalized so 
that each quark or anti-quark field contributes 1/N C to the quantity kjjy, when J is the baryon 
current of (2A). If we equate the values of the left side of ( |2.24 ) at scale /i and in the UV limit, 



we find 

Ka) = buv + i(g)~kuv - (l - - 7(2)) [%) - kuv] • (2-25) 

This result can be checked against the explicit 2-loop calculation M of the flavor central 
function b{g) and it agrees. We do not discuss this here, because we are more interested in non- 
perturbative application to electric SUSY QCD in the conformal window, 3N c /2 < Nf < 3N C , 
where there is evidence that the theory flows to a non-trivial fixed point g* Q. Although g* can 
only be calculated at the weakly coupled end of the conformal window (where Nf = 3iV c (l — e)) 
we know that the /3- function fl2.5|) vanishes, and this is enough to give the exact infrared limit 
of the anomalous dimension 

1IR = (l ~ (2.26) 



N 



f 



Then fl2~25p becomes 



bm - b uv = imkuv, (2.27) 

which is our first non-perturbative result for the flow of a flavor central charge. Since gauginos 
do not contribute to flavor current correlators and the quark and anti-quark contributions to 
the combination R + K/3 cancel, it follows that buv = —kuv f° r correlators of all of the flavor 
currents in (p.l|). For the baryon current, buv = 2Nf/N c and we have the total flow 

bm -b UV = eU- , (2.28) 

which is positive in the entire conformal window, contrary to c-theorem intuition. 

We now turn our attention to magnetic supersymmetric QCD with gauge group SU(N' C ), 
with N' c = Nf — N c . The matter content consists of Nf flavors of magnetic quark and anti-quark 
superfields qf(x,9), q^(x,#) plus the gauge neutral Nf x Nf meson superfield Mj. There are 
conventional gauge interactions and a cubic superpotential W = fqMq. The dual theory has 
the same flavor group SU(Nf) q x SU(Nf)~ x U(1)b, but for simplicity we shall consider only 
the baryon current n 

11- 

JpL = N 1 2 (^7^75^ - ^qlfxl^q) ■ (2.29) 

We also need the separate Konishi currents of the quarks and the mesons 

4 9) = \ (Vws^ + i><n^i>o) , (2.30) 



For simplicity we do not write the scalar contributions to the currents in eqs.(2.29-2.31) 
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while the fermion content of the R current is 

R, = ^ 7/i75 A a - \{K^ + KM). (2.31) 

The divergences d^K^ and d^K^ have classical contributions, since the superpotential is 
not invariant under the relevant U(l) transformations, and duKjfi also has an anomaly. The 
internal i?-current anomaly also involves FF and there is a superpotential contribution (II], 113] . 
It is significant for our analysis that there is a unique combination |l0| of these currents which 
is classically conserved and anomaly-free, namely 

S, = R, + ~ (l - ^ - 7ff ) (4* " 2K l M) ) ~ I ^ + 1") 4 M) ( 2 - 32 ) 



The combination — 2K^ is conserved classically. Its coefficient in ( 2.32j ) is again the 



numerator of the NSVZ gauge (3- function ( |2.5| ), and the coefficient of K^" 1 is essentially the 
Yukawa /3-function j3f = f(2 r y q + 7m)- 

The previous analysis must be generalized to include a superpotential. There are only trivial 
changes in the derivation of the relation between the external trace anomaly b and the central 
function b(g(l/x), f(l/x)), and the result that they coincide at fixed points is unchanged. The 
analysis of the flow of flavor central functions is also easily repeated. The third equation of 
( p. 23 ) is replaced by the pair 



and one finds 



(d^Kft) = - JL ~ k M {g{ ^ + ... (2.33) 

<<yrw> = - JLpo GKm), /(m)) b, v b, u + ... 



b = b uv + jfijl + iMk^ (2.34) 

2A-( M ) _ k {q) + 2k {MY 

Iv/jy -|- ^l^jjy 




+ {2 lq + 1M )[k^-k < : M) 



uv 



in which the subscript UV denotes a quantity evaluated in the ultraviolet from 1-loop graphs, 
and other quantities are evaluated at scale \i. One can also check, as in the electric case, that 
this formula agrees with the two-loop calculation. 

We now go to the IR fixed point where the coefficients of the last two terms of ( p. 34 ) vanish, 



since they are related to f3 g and (3f, respectively. This gives the IR values of the anomalous 
dimensions 



10 



and we obtain from (2.34) 



(7^ + 7^), (2.36) 



which is a non-perturbative formula for the flow of flavor charges in the magnetic theory. For 

r $l = 2N C /N' C and fcg# 



the baryon current, buv = ~^uv = ^N c /N' c and kyy = 0, so the flow of the baryon central 



N f \ 2N f -3N C . 

'"' '•^-WI) =2 ^nT (2 ' 37) 



charge is 



which is again positive throughout the conformal window. From ( |2.28| ) and ( |2.37j ) one can 
check the equality of the IR values of the central function, namely bm = 6. This is not a 
new confirmation of duality, since it is the same 't Hooft anomaly matched in electric-magnetic 
duality gj. 

We briefly discuss the question of constancy of the flavor central charge on conformal fixed 



lines [18 1, confining our attention to the N = A and N = 2 cases. The matter content of the 
N = 4 theory is three adjoint chiral multiplets, while for N = 2, with G = SU(N C ), Nf = 2N C 
fundamental hypermultiplets are required to produce a fixed line. In each case there is a 
particular cubic superpotential, and the condition for the fixed line is a linear relation between 
g and the Yukawa coupling /. It is straightforward to generalize the previous derivations 



and show that formulas of the type (2.27) and ( |2.36| ) hold for the N = 4 and N = 2 cases, 



respectively, with [TV-subscripted quantities calculated from the relevant 1-loop graphs. To 
establish marginal constancy at the non-perturbative level, we need only observe [18] that 



anomalous dimensions vanish on the fixed line, so we have b =constant. At first thought one 
might expect that the value of the central charge depends on the single coupling, say g, that 
remains on the fixed line; however we see that b is determined by the free field content of the 
theory. 

An even simpler argument is to note [jl0| that Sn and coincide at infrared fixed points or 
fixed lines. But is anomaly-free, so we have {d fl S^J v J p )uv = {dfiS ll J l/ J p )iR = {d^R^JyJp) ir. 
Hence, the central charge b on an infrared fixed line is equal to the anomaly coefficient — sq and 
thus depends only on the free-field content of the theory and not on any coupling constant. 



3 An alternative approach 

In this section we compute the flavor central charges using a somewhat different approach. 
This method uses the correlation function b(g(l/x)) of ( p.6[ ) in a more direct way, and does not 
require an explicit regularization of the singularity at x = 0. This simplifies the analysis. The 
method is quite general, but we work with the baryon number current of the electric SUSY 



QCD, given in (2.1) for simplicity. 

Let us consider the correlator (|2.6|) again and compute the scale derivative \idjd\i{ J^(x) J v (0)) 
at fixed bare coupling go = <?(//), which gives by definition the correlator (J^x) J v (0) J d 4 z @(z)). 
In the presence of an external field coupled to J p and in quadratic approximation the matrix 
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element < / d z @(z) > has the form 

< J d'z e(z) > = / (0i /GK*0) • \B^{k)B, v (-k) = (3.i) 

d A k 1 „„„,.,„ , , , / „ 1 /" ./i 1 9 

1 „„„,.,„ , f „ I f .A d 

W) 



l -B^{k)B, v {-k) (fo-^J d^^^b{g{l/x)) 



The above integral is convergent at x = because of asymptotic freedom and the RG equations. 
In eq. Q3.1| ) we explicitly introduced a constant fo associated with the regularization of the 
correlator (|2.6|) at x = 0. This contact term is irrelevant for the present analysis. However it 



can be easily determined. Indeed, in the UV limit, i.e. \k\ — > oo, the x-integral in eq. ( 3.1 
vanishes, and hence fo = fuv = /(<?(&)) ||fc|-»oo • ^ n particular in the case of the baryonic 
current / = Nf/4ir 2 N c . 

The IR limit corresponds to \k\ — > 0, i.e. fjR = f{g{k)) |ifc|_>o- Thus from (|3.1[) one gets 

1 f°° d 1 

fuv - fm = -5-2 / dlogx 2 ^b(g{l/x)) = -^(b uv - b IR ). (3.2) 

8tt z Jo ologx z 8ir z 

We turn now to the computation of the function f(g(k)). The simplest way is to observe that 
the B^J^ coupling is essentially that of SQED, with external "photon" field B^ and coupling 
constant e. In SQED the conformal anomaly can be computed as a scale derivative fid/dfj, (at 
fixed values of g(n), e(/i)) of the effective action for the external field Bn 

™ = Jl&4*® B ~ ik) "~ i - k)t (3 - 3) 

where e(k) is the effective running coupling constant. We have 



< 9 W > = rg-m - - J B( n( -k) ■ 



dp J (2tt) 4 v ' ^ v y 2e 3 (fc) 



A? 7 (o.^^r^^-^rll-T^))]. (3-4) 

This expression follows from the general result for the N=l supersymmetric QED /3-function 
Pe = n (1 — 7)e 3 /87r 2 given in ref. [11], 7 being the anomalous dimension of the matter chiral 
superfields and n standing for the number of flavors. In the present context we have n = N c Nf, 
and e — > because the QED gauge field B^ is external. Therefore, 

/(#)) = Jk [1_7(#))1 ' (3 - 5) 

gives an exact formula for f(g(k)) as a linear function of the anomalous dimension. (The 
relation ( |3.5| ) is not scheme independent, but rather holds in the scheme of |l(J.) The inverse 
Fourier transform of ( |Q| ) then gives fi-^ (b(g(l/x))/x^) , and b(g(l/x)) can be easily computed. 
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Since r y(g(k)) |fc^oo= and r y{g(k)) |fc^o= 1 — 3iV c /iV/, from eq. (3J3) immediately follows 
that 



bm - b uv = 87r 2 (f IR - f uv ) = 2 • ^ • ( ^ - 1 1 . (3.6) 



We therefore recover the result of Section 2. Note that we did not use the current of ( [2.4D . 

It is worth observing that the SQED /3-function is also implicitly present in the treatment 
of Section 2. To show this we consider a theory with gauge group extended to SU(N C ) x U(l) 
with U(l) coupling constant e. This theory has the operator trace anomaly 

° " 32^ (W + lQ^Nl {Bfiv) ' (3 - 7) 

where the coefficient of (B^) 2 is (3 e /2e 3 N 2 and 7 = j(e(fi), g(fi)). In the limit e-»0we can 



compare (|3.5|) and (^J) with (|2.9| ) and identify 



2/3 e ^(1-7(0,^)) 

^ ^ivf = i^jv? = f ^ )] (3 ' 



It is instructive to compare {3A) and (2.21). Consider first (2.21). The correlator ( J M J u J O) 
is renormalization group invariant. In momentum space the r.h.s of (2.21) may depend only 
on the running coupling g(k), but not on g(n). Therefore, the //-dependence of q is compen- 
sated by that of the correlator {J^Jy J{F® U ) 2 } which has a perturbative expansion starting with 
g A (/j,) log k 2 /fi 2 . On the other hand the form of ( |3.1| ) is explicitly renormalization group invari- 
ant. The local part q of (2.21) was identified above as q(g{fi)) = f(g(fi)). The non-local part of 
([O]) corresponding to the contribution of {J^Jy J(F^ U ) 2 } in (2.21), is then 

/($(*)) " /G?(M)) ~ [3iV c -N f (l- 7(0(M)))] • [g\ri log k 2 /^ 2 + ...]. (3.9) 
This is proportional to (3(g(fi)) because of the RG equation 

a* l^(iAO) = pw)) Qg^jr^M) (3-io) 

which can be used in ( |3.1|) , and perturbation theory gives the leading power g A (n). In this way 
we establish a correspondence between the approaches of Section 2 and the present section. 

Another method of calculation is to consider first the anomalous superspace operator equa- 
tions in the absence of the external field (see, for example JO], |H|) 

D^-^(l-^-^D a W 2 = 0, D 2 K-^W 2 = 0, (3.11) 

where J a a is the supercurrent, K is the Konishi operator. By combining eqs. ( |3.11| ) one can 

get 
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which is the superspace version of (|2.4j). In the presence of the external field B^ the operator 
D a J a a — (l — 7§y — 7(5)) D a D 2 K does not vanish but is proportional to the unit opera- 
tor times a local functional of the external field. Hence, for the component matrix element 
containing the stress tensor we have (0 oc [D a , D a ] J a a\e=o, ^matter ot {D 2 , D 2 ]K}e 



te=o 



< 



d 4 z 



Q(z) 



1 



N 



f 



7(5) L matter (z) 



> 



-k). (3.12) 



The constant no does not depend on the external field. It is also independent of the renor- 
malization group scale \i since the operator on the left hand side is renormalization group 
invariant. 

(It is worth noting that the above statement is a consequence of a general result. We refer 
to an operator O that has vanishing matrix elements between any Fock physical states as a null 
operator. This is described by the operator equation O = 0. Thus, D a J a a — ^(1 — 3N C /Nj — 
r y)D a D 2 K is a null operator. In the presence of an external field a null operator may not 
vanish. In general the r.h.s. of the operator equation may be a linear combination of various 
local operators Oi with the coefficients yi being non-trivial local functionals of the external 
field, i.e. O = J2i V%Oi- In the simplest case the O operator mixes only with the unit operator, 
i.e. the operator equation reads O = y ■ 1. In such a case the matrix element of O is just a local 
functional y of the external field. Note that the one-loop form of the 't Hooft external anomaly 
for a current is a particular consequence of this fact. Indeed, assuming that does not 
have any internal anomaly its anomalous dimension is vanishing. Therefore the matrix element 
of d^J^ is renormalization group invariant, i.e. it may depend only on the running coupling 
g(k), where A: is a momentum of the external field. On the other hand, as mentioned above, y 
is a local functional of the external field, and hence it does not depend on g(k). Therefore the 
matrix element in question is exactly one-loop.) 

For technical reasons it is convenient to consider a particular momentum mode k and define 
the reduced matrix elements <C • • • 3> obtained by dividing the conventional matrix elements 
(• • •) by (1 / ^jB^ (k)B ^{—k) . The reduced matrix element of the operator J d 4 z [@(z) — (1 — 
3N c /Nf— r y(g)) L ma u er (z)\ (which is just equal to uo) does not change under the renormalization 
group flow with respect to k 2 . Therefore we have 



A/ = f(g(k)) || fcH0 - f(g(k)) |, fcHoo =« / <fz (1 - 3N C /N f - 7(g(fi)))L matter (z) »|| fcH „ 



< / (fz{\ - W c /N f - j(g(n)))L mat ter(z) M\k\ 



(3.13) 



The matrix element <C / d 4 z (1 — 3N C /Nf — 7(5)) L mat ter(z) ^> is proportional to 1 — 3N C /Nf — 
r y(g(k)) pf and hence vanishes in the infrared, i.e. at k 2 — > 0. Thus we have, precisely as before, 



< / d A z (1 - W c /N f - 7(5)) L matter {z) >|| fc |. 



Nf 




(3.14) 
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where the result follows from the identification of L ma tter with the 6 2 6 2 component of the 
Konishi superfield, and the (one-loop) Konishi anomaly. 

The story becomes more complicated for the gravitational central charges which are con- 
sidered in section 4. In this case one may try to analyse the matrix element < J d 4 z Q(z) > of 
the stress tensor in the presence of an external gauge field field coupled to the current in 
order to compute the flow, say, of the c charge. However, the coupling of the current to the 
"QED" vector field breaks supersymmetry and this "QED" /3-function is not simply related 
to the anomalous dimensions 7. A direct way to compute the flows of the c and a charges is to 
consider the matrix element of the divergence of the R^ current in the presence of an external 
gauge field coupled to the i? M current. In this way it is easy to rephrase the considerations of 
Section 4 in the formalism of this section. 

4 Gravitational central functions 

As discussed in || the correlator of two stress tensors can be written in the form 

, m , , m 1 „ c(q(l/x)) f(\nxu,q(llx)) . 
(T^(x)T p<7 (0)) = -— V V J + x l > 

where tt^ = (d^d v -5^ v u) and n^ pcr = 27^7^-3(7^7^+7^0-7^) is the transverse traceless 
spin 2 projector. One of the objects of primary concern in this section is the flow of the central 
function c(g(l/x)) that we will relate to the coefficient of the square of the Weyl tensor in the 
external gravitational trace anomaly. The second object of interest is the coefficient a of the 
Euler term in the trace anomaly. We will give non-perturbative formulas for the flows of these 
quantities. The internal trace anomaly is responsible for the second tensor structure in (|4.1| ). 
It is proportional to /3(g(p)), and thus vanishes at critical points. 

Since T fiu (x) and R^(x) are both components of the supercurrent superfield J a a(z), z = 
(x,8,9), there is a relation between the (TT) and {RR} correlators. For a critical supersym- 
metric theory, they are both contained in the supercorrelator 

(Jaa(z)J p$ (0)) ~ C (4.2) 

given in Section 5 of ||] where the notation is explained. This means that the central charges of 
the TT and RR OPE's are given by the same constant c which is a fixed point value of c(g(fj,)). 
Off criticality there is a more complicated relation between (TT) and (RR) that is not required 
for our work. 

In Section 2 we introduced a source for the flavor current J^{x) in order to show that 
the coefficient of the external trace anomaly coincides at fixed points with the central function 
b(g(l/x)) and related the trace and d^R^ anomalies using global supersymmetry. In the present 
case sources for T^ u and R^ take us out of the realm of global supersymmetry and require 
external field supergravity where the anomaly situation is more complex, as we now discuss. 
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We introduce the background metric gu V (%) and source V^x) for the i?-current. In these 
background fields the trace anomaly of a critical supersymmetric theory has the form 

6 = 16^ {W ^ ? ~ 16^ + ^2 K ■ ( 4 - 3 ) 

Here W^ vpa is the Weyl tensor and R pvpa is the dual of the curvature tensor, the second term 
being the Euler density; is the field strength of V^. The coefficients of the iW^pa) 2 and 
iV^y) 2 terms are related while that of the Euler density is an independent constant. In a free 
supersymmetric gauge theory with N v gauge and N x chiral multiplets, the constants are [|l^] 



n\ - (3iY„ + N x ) , a uv = ^ (9N V + N x ) . (1.1 ) 



Off criticality there are additional terms in O that are proportional to and do not 

contribute to the total flow, and the central charges depend on the coupling, i.e. c = c(g(n)) 
and a = a(g(/j,)). 

In the superspace description of these anomalies, the external metric (actually the vierbein 
e^) and current source (the supergravity axial vector auxiliary field) are contained in a single 
superfield H a (x,6,9), and the trace anomaly and the dpRp anomaly are components of the 
superfield equation 

D«Jaa = D a J (4.5) 

where the chiral superfield J (the supertrace) has the form 

J=^(cW 2 -aZ c ) (4.6) 

in terms of the superWeyl tensor W 0j g 7 and the chirally projected superEuler density (see 
Appendix A for details). 

We need the d^R^ anomaly. Although it is fairly straightforward to compute it as the theta 



component of (4.5), and we will do this in Appendix A, we present here an indirect argument 
that uses the general structure of (|4.6|) but does not require any superspace technology. The 
d^R^ anomaly will have the form d^R^, = (uc + va)R a b c dR abcd + (wc + za)VV with model 
independent coefficients u, v, w and z. It is thus sufficient to compute in a free supersymmetric 
gauge theory in order to find these coefficients. Because of the ratio —1/3 between the i?-charges 
of gauginos and matter fermions, we know that uc+va can be obtained from the (TTR) triangle 
graph as a pure numerical (i.e. 1/247T 2 ) multiple of 3^ — N x , and wc + za can be obtained 
from the (RRR) triangle as a numerical multiple of 27 N v — N x . Using the values for c\jy and 
ajjv in ( fOl ) we find 

^(vW) = ^iW^ + (4-7) 

The next step is to note that the W 2 anomaly coefficient of a non-critical interacting theory, 
c(g(fj,)), can be related to the central function c(g(l/x)) by the same argument as in Section 2, 
which in fact is exactly the argument given in ||. (The Euler term of gives no contribution 
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to the integrated trace anomaly.) The functions c(g(/i)) and c{g{(j)) coincide at fixed points of 
the flow. 

We now proceed to the calculations of c and a in the electric SU(N C ) SQCD, where N v = 
- 1 and N x = 2N c N f . To achieve this goal, we can apply arguments of Section 2 to calculate 
the flow of the combinations of anomaly coefficients c(g(fi)) — a(g(fj,)) and 5a(g(fi)) — 3c(g(fj,)) 
from the three-point correlators (TTR) and (RRR), respectively. 

The anomalies of the correlators {TTS), (TTR), and (TTK) are summarized by the equa- 
tions 



{d^W)) = JL (c(g(v)) - a{g(n))) e^R^ Xr R%, 



(d^K*)) = -^k(g^))e^R pvXT R%, (4.8) 

where we omitted the non-local 0(f3(g(fi))) terms. In the ultraviolet limit the quantity k(g(fi)) 
can be obtained from the 1-loop triangle graph as k — ► kuv = —NfN c /S. Similarly c(g(/i)) — 
a(g(fj,)) — > cuv — o-uv = ~Jq(^c — 1 — f-^/^c)- It is now immediate to write the analogue of 
( P^D which is 

1 if 3N C \ 

c(g) - a(g) = c uv - a uv + -j(g)kuv ~ 3 ( 1 ~ ~J^J ~ I (%) ~ k uv)- (4.9) 

To discuss the (RRR) anomaly, we observe that the amplitude of the triangle graph for 
the contribution of one Majorana spinor with current J M = ^ 1^7^75 ^ has the Bose-symmetric 
anomaly 

d 1 d d 16 

-Q^(Jfi{ x )Ju{y)J P {z))uv = ~J^ e » v P°far&jj~ s ( x ~ z ) 6 ( y - z ) = -Q\u(x,y,z) (4.10) 

(see [21 1 for a recent discussion of the anomaly in x-space). We then write i? M = 5 M + ^(7 — 
jiji)K p where Su is the internal anomaly- free current ([O]). The various Bose symmetric 
contributions to (RRR) then have anomalous divergences that can be written as: 

-^-(R p R v R p ) = [5a( 5 0u)) - 3c(< ? 0tx)].V 

\(S p S v K p ) + (S^KvSp) + (K p S u S p )] = Zk^gin^A^ 
^- \(S p K v S p ) + (K p K u S p ) + (K p S u K p )) = 3k 2 (g(^)A tlv 



d 



(K p K v K p ) = MsG"))-V (4.11) 
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where the tensor indices n, u, p, are associated with coordinates x, y, z, respectively; we again 
omitted 0((3) irrelevant nonlocal terms. The anomaly coefficient S3 is scale independent while 
other anomaly coefficients depend on g(fi). 

Proceeding as before and evaluating S3 in the UV limit from triangle graphs, we obtain 

5a — 3c = 5auv — 3cuv + h, (4-12) 

where 

h = jmh uv - lm\k 2 uv + im^huv + (7 - im)fa 

+{l-liR) 2 \k 2 + (n-iiRf^h. (4.13) 

The subscript UV indicates the free field value, while the other quantities are evaluated at the 
scale fi. We have 

1 9 9 N c 9 / N c \ 2 , 

kuv = - V6 N x , k 3UV = -N x , k 2UV = --N x -, huv = T - Q N x ^-j . (4.14) 

Finally, collecting (g^) and ( jig) , we find the non-perturbative results 

5 5 1 

c = c uv + 77(7 - lm)k + -jmkuv + -zh, 
b o 2 

a = a uv + -(7 - jm)k + -imkuv + zh. (4.15) 

Let us compare our results with two-loop calculations. The work of || shows that a (g) has 
no 2-loop corrections and in JJJ the result 

1 < P(g)\ 



c(g) = - [3N V + iV 7 - 7I + N v ^j (4.16) 

for c was obtained. We see that the two-loop correction is the sum of a contribution proportional 
to 7 and a contribution proportional to (3(g) /g. Formulas ( [4.15 ) show that the coefficients of 
the 7 terms are 

-kuv + -huv ~ ^Jmk2uv + jglmfauv = --^N c N f , 

\ k uv + 2 fa uv - ^Jmk2 uv + jgJiRfa uv = 0, (4.17) 

for c and a, respectively, in agreement with the two-loop results of || and ( 4.16j ). Note that 
7^ = N c Nf^f in the present notation and 7 is the anomalous dimension of the gauge invariant 
QQ which is twice the value of the field anomalous dimension used in [J5J . The comparison of 
the coefficient of the (3(g) jg terms, instead, is more subtle and requires the precise knowledge 
of the two- loop corrections to the functions k, ki, k% and k^. For example, a two-loop correction 
to k is expected from the photon chiral anomaly in an external gravitational field |20||. 
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In the infrared limit 7 — ji R vanishes (it is the numerator of (3{g)) and the central charge 
flows are 



N c N f ( N c Ni \ N c N f 9 N c \ , 

cm - c uv = -^- 1IR (3— + 9^ - 4j , a IR - a uv = —^f IR (2 + 3— J . (4.18) 

We can check agreement with the two loop results also from these formulas, but only in 
the weakly coupled region ^i R <C 1. Again, the coefficients of ji R are the expected ones. The 
fact that a is two-loop uncorrected is exhibited by the appearance of jj R in the expression of 
a-iR — auv- The term proportional to f3(g)/g in ( 4.16|) , on the other hand, cannot be reproduced 



in the flows ( [4.18 ): although it is nonvanishing at the order g 2 , it is cancelled at criticality by 



the higher loop corrections. 



In the considerations above we have discussed the central functions along the RG flows. If 
one is only interested in the fixed point values of c and a , one can use a shortcut because the 
coefficient (1 — 3N C /Nj — 7) vanishes at the IR fixed point, so Rp and effectively coincide 
there. It is again crucial that the internal anomaly-free current Sp of (|2.4| ) has one-loop exact 
external anomalies, so that the IR values of the it^-anomalies coincide with the UV values of 
the corresponding ^-anomalies, which makes them computable. Thus, for example, 

d d 

{Rfi{x)R v (y)R p {z))iR = - — {S ll {x)S v {y)S p {z)) IR 



() (S^x)S u (y)S p (z)) uv (4.19) 



dxp 



and the last anomaly can be obtained from the one- loop graphs of S^ v = Rp+\{l-3N c /N } )Kp 



with Rp and Kp defined in ($L2\). We find 

5a IR - 3c IR = ±(n*-1- 2N c N f ^) ) (4.20) 

In the UV limit, ■^-(R(x)pR(y) u R(z) p ) is determined by the one-loop fermion triangle graphs 
of Rp in (O) which yield 



9 / 9 2N c N f \ 
ha uv - 3c uv = — (^ c 2 - 1 - • (4.21) 

A similar argument applied to (TTR) can be used to determine the flow of c — a. One finds 

ci R -a IR = —(N* + 1) 
lb 

cuv -a uv = -^ U 2 - 1 - ~N f N c ) (4.22) 
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(The triangle graph contributions in (4.20-42^) are normalized consistent with (4.4).) The 
flows of 5a — 3c and c — a agree with the infrared limits of ( |4.12| ) and ( [4.9D , respectively. 

To conclude, the IR values of the two gravitational central charges in the electric theory are 



cm 



7Nt 




am 



2Ni 



1 




(4.23) 



These results will be discussed in the next section. 



5 Discussion 

The techniques developed here are quite generally applicable to any supersymmetric gauge 
theory, provided the theory flows to an IR fixed point and the gauge /3-function has the general 
structure of the NSVZ formula, so that anomalous dimensions can be determined at the fixed 
point. Specific applications have been made to the electric N=l SU (N c ) series in the conformal 
window 3N c /2 < Nf < 3N C and their magnetic duals. 

It should be noted that our formulas for central charge flows do not give new tests of 
duality because the current whose anomalies agree in the electric and magnetic theories 0] 
coincides at the IR fixed point with the Rn current for which the anomalies are various linear 
combinations of the IR central charges. We derive here the flows cm — cjjv and am — a\jy in 
the magnetic case by this shortcut. The IR values are just those given in ( 4.23] ) for the electric 



theory, and the UV values are determined by the free-field content. In this way we find, in the 
magnetic theory, 

°« " <*" -k (> " \w) ( 9 f " SN < + SN f + H • 

a '*-WV = ~\l-~) (3JV c 2 + 4JV c JV / + 3JVf). (5.1) 

In these formulas, 7^ is an overall factor and, again, the correct two-loop results in the weakly 
coupled region of the magnetic theory, 7^ <C 1, are reproduced. Note once again the appear- 
ance of (7g R ) 2 in ajR — ajjv- 

Let us discuss our results for central charge flows from the viewpoint of c-theorem expec- 



tations. The results ( |2.28| ) and (|2.36| ) for central charges of flavor currents in the electric and 
magnetic theories both satisfy bm — bjjv > in the entire conformal window, so there is no 
c-theorem for flavor. 

For the central charge c of the TT OPE and the Euler anomaly coefficient a, there is 
earlier work by Bastianelli p2| for SU(N C ) SQCD with Nf = and in the confinement range 
N c < Nf < 3N c /2. The IR limit is a free theory of massless excitations of independent gauge 
invariant composite operators. In an empirical approach, Bastianelli computed the free field 
values in the IR and UV and found both cm — c\jy < and am — a\jy < 0, thus suggesting that 
a c-theorem holds in supersymmetry for the central charges c and a. Our techniques permit 
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the extension of these c-theorem tests to the conformal window where there is an interacting 
IR fixed point. 

The results for c show both positive and negative flows. In the electric theory cm — cjjv 
is negative near the lower edge of the conformal window, but positive near the upper edge. 
In the magnetic theory cm — c\jy > in the entire range. We conclude that there is no 
supersymmetric c-theorem for c. Discussions about this issue have already appeared in the 
literature in the domain of weakly coupled theories, both supersymmetric , for example in [fjj 
using formula ( 4.16| ), and nonsupersymmetric |23}| . 



The story of the Euler central charge a is very different, since we find am — ajjv < for both 
electric and magnetic theories in the full range of the conformal window. We believe that this 
is strong (because nonperturbative) new evidence in support of an a-theorem (i.e. a c-theorem 
for the Euler anomaly) and thus of irreversibility of the RG flow in quantum field theory. This 
is the first time that these ideas have been tested in strongly coupled theories. 

One can also consider the central charge flow when quarks masses are introduced or the 
Higgs mechanism is used to give masses to some gauge fields and their superpartners. In either 
case, the IR values of the central charges change, while the UV values are the same as before. 
One finds, from ( [4.23 ) and the free field UV values, 



am ~auv = ^[ 18iV c 2 - 27^ - 9N' C 2 - 2N' c N' f ) . (5.2) 



Here the IR values N c and Nf are assumed to be in the conformal window, in order to have 
an IR fixed point, while the UV values N' c and Nf are subject only to the conditions N' c > iV c , 
Nj > Nf and Nj < 3A^. The condition am — auv < is satisfied in this more general situation, 
where some fields are integrated out in the flow. 

Let us discuss other positivity conditions suggested by the simple physical intuition that 
the central charges count degrees of freedom, variously weighted according to the external field 
one is using (gravitational, flavor, etc.). In particular, one expects the following inequalities to 
hold, 

C ^- ' a ^- ' 9iv>- ' 5iv7- ' aivT- ' 9aT-°- (5 ' 3) 

The first inequality is rigorous, since c appears in the < TT >-correlator. The second inequality 
states that a also should count degrees of freedom and has been recently proposed in j^4j. The 
other inequalities express the expectation that adding matter fields (by varying Nf) or enlarging 
the color group (by varying A^ c ) should increase the number of degrees of freedom. 

The inequalities ( |5.3| ) obviously hold in the UV and we would like to discuss them in the 
IR. Imposing them on our results flQjp we find the following restrictions, in the large N c and 
Nf limits, 



N}>-Nl N 2 f >-Nl n.r., N 2 f > -iV c 2 , n.r, N 2 f > 3iV c 2 , (5.4) 
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where n.r. means "no restriction". The first inequality, as we said, is rigorous and so puts a 
limit on the region where the IR fixed point exists. This region contains the full conformal 
window, in agreement with electric-magnetic duality. We see that below a certain value of Nf 
our treatment necessarily breaks down, which means that no IR fixed point exists. This is 
expected, since, for example, pure supersymmetric Yang-Mills theory (Nf = 0) has no IR fixed 
point, according to the NSVZ exact /3-function. 

The other inequalities of (|5.3| ) have not been proved and so our discussion is purely specu- 
lative. One can observe that all but the last one are satisfied in the entire conformal window. 
So, if one assumes electric-magnetic duality, then the physical intuition that we suggested is 
not completely correct and there is a region of the conformal window where enlarging the gauge 
group decreases the value of a in the IR. We do not have a way to resolve this puzzle at the 
moment, but we believe that these remarks are relevant in order to understand the nature of 
central charges better. 

A final comment concerns the relation between central charges and 't Hooft anomalies. The 
't Hooft anomalies are quantities that are constant at all energies, so computable at the free 
fixed point. They can be regarded as the invariants (or "indices") of a quantum field theory. 
The central charges, on the other hand, have a direct physical meaning and their RG flow is 
nontrivial. In supersymmetric theories the RG interpolation problem can be solved because at 
the IR fixed point the central charges are related to the 't Hooft anomalies. We now show that 
all the known 't Hooft anomalies can be expressed in terms of the central charges. 

The 't Hooft anomalies that contain at least one vertex U(l)s are related by supersymmetry 
to appropriate terms in the trace anomaly and therefore to primary central charges. We know 
that the U(l)s anomaly is proportional to c — o, the U(l)g anomaly is proportional to 5a — 3c 
while the U(\)sGG' anomaly is related to the flavor central charge b, G and G' denoting flavor 
groups other than U(l)s- Other G\GiG 't Hooft anomalies do not contain the vertex U(l)s- 
They are related to certain secondary central charges c'. Using the construction of Section 3 
of ||, let us consider the G-channel of the {G\(x)G2(y)Gz[z)G±{z)) four-point function, where 
the limit \x — y\, \z — w\ <C \x — z\ is taken. These secondary central charges are rather special, 
since the intermediate channel is again a conserved current. The value of c' is simply expressed 
in terms of a primary central function &g,g and two 't Hooft anomalies. More precisely, 

<&<;G = (GiG 2 G) b GiG (G, G 3 , G 4 ) (5.5) 

(GiGjG) denoting the value of the appropriate triangle anomaly. Therefore the set of primary 
and (special) secondary central charges contains the set of 't Hooft anomalies. An implication 
of this remark is that the central charges in question coincide in both the electric and magnetic 
theories. The secondary central charges that are not included in this set, on the other hand, do 
not correspond to known 't Hooft anomalies, but they should also coincide in the electric and 
the magnetic theory. One possible extension of the analysis of this paper is the study of these 
new central charges, with the purpose of testing non-abelian electric-magnetic duality. 

It would be interesting to extend our investigation to nonsupersymmetric theories. The 
properties of supersymmetry have been intrinsically used in our derivation. In particular, 
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crucial roles were played by the exact expressions for /3-functions and the relation between the 
trace anomaly and the chiral anomaly of the i?-current. In view of this, the nonsupersymmetric 
generalization is nontrivial. A relatively simple case to start with is pure Yang-Mills theory. 
The central functions c(a) and a(a) should flow to zero in the IR, since the theory is expected 
to have a mass gap. 

Finally, the other interesting open problem is to rigorously derive the a-theorem (at least 
in the class of supersymmetric theories) and definitively prove that the RG flow is irreversible 
in quantum field theory. 
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A Appendix 

We present in this Appendix a superspace derivation of the coefficients c — a and 5a — 3c 
appearing in ( [4.7| ) . 

For any matter system coupled to (background) supergravity the supercurrent and super- 
trace are defined by 

ST 7 _iL (A u 



where T is the (classical or quantum) action, H aa is the supergravity vector prepotential super- 
field, and (j> is the superconformal (chiral superfield) compensator. The superfield H aa contains 
the vierbein and, as its last, 6 2 8 2 component, the axial vector auxiliary field denoted here by 
V aa . The supercurrent and supertrace are related by the (local supersymmetry) conservation 
equation 

V d J aa = V Q J (A.2) 

If the theory is superconformal the action is independent of the compensator and the supertrace 
vanishes. Otherwise, J contains the trace anomaly, the supersymmetry current 7-trace anomaly, 
and the chiral anomaly for the i?-current. The i?-current appears as the first component of the 
supercurrent or, equivalently, as the term in the action that couples to V aa . 

For supersymmetric matter systems (scalar or vector multiplets) the supertrace has been 
computed in terms of so-called super-64 coefficients in refs. [po] , |26| ( however, we use in this 
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Appendix the conventions of Superspace |27j|). One finds 
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where the numerical coefficients c, a take the values 

Scalar multiplet : c 
Vector multiplet : c 

Here W 2 = \W aPl W a ^ is the square of the superWeyl tensor, while 

E c = W 2 + (V 2 + R){G 2 + 2RR) 
is the chirally projected superEuler density. More precisely 
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give the Euler number and Pontrjagin number, respectively. The chiral superfields W a p^ and 
R, and the real superfield G a a are the three superspace curvatures. In our conventions G 2 = 



1 /~iaa/~< 



We are interested in the component trace and i?-current anomalies for the matter system 
in a background gravity and axial-vector auxiliary field i?-current source. The trace of the 
stress-tensor is obtained from the 9 a 9 a component of the supercurrent: 



9 = -[V a ,V Q ]J Qci | e= o = -(V z J + V z J)\ e=0 (A.7) 
while the divergence of the i?-current is obtained from the first component of the supercurrent: 



1 



iV a R a = ^{V a ,V a }J a a\e=o = (VV- V 2 J)\ 0=O 



One can obtain the corresponding component anomalies by brute force ^-expansion or, as we 
shall do here, by definining components by projection and exploiting the Bianchi identities of 
the theory. In doing so, we will obtain component results expressed in terms of the component 
curvature tensor R a bcd and the axial vector field strength, V a b- However, one further step is 
necessary: the standard supergravity constraints [27] lead to component covariant derivatives 
with nonzero torsion, due to both the gravitino fields (that we set to zero here), and to the 
axial vector auxiliary field. Before reading off the component anomalies, we must separate the 
curvature in terms of the ordinary Einstein curvature and additional terms, again proportional 
to the field strength V^. 

We obtain the component results as follows: the W 2 term in the supertrace leads to 



^V 2 W 2 ± ^V 2 W 2 = jV s W a ^V s W aPl - jW a ^V 2 W a ^ ± h.c. 



(A. 
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evaluated at 9 = 0, the db sign corresponding to either the trace or the chiral anomaly. Since 
we are restricting ourselves to a bosonic background the second term (proportional at 9 = to 
the gravitino field strength) may be dropped. We write each factor in the first term as the sum 
of symmetrized and antisymmetrized (in 5 and the already symmetrized a/37) indices, i.e. 

V 5 WVy = ~ [V {s W al3l) + C a5 V X W XM + C ps V x W Xia + C 7(5 V A W Xaf3 ] (A.10) 

For the last three terms we use the Bianchi identity V A W\/3 7 = f^^G ^. Evaluating at 
9 = 0, the (totally symmetrized) first term gives the self-dual part (in spinor notation) of the 
component Weyl tensor, = ^jV '(fW^i (but we emphasize again that it still contains 

some torsion pieces). The second term, using the definition of the axial vector auxiliary field 
Va = C a |6i = o leads to the self-dual part of its field strength, Vg 7 = \d,p^vP . 

For the G 2 terms in the superEuler density, dropping irrelevant terms, we have 

(V 2 V 2 ± V 2 V 2 )G 2 = --[V a V d G^V a V d G^ ± V 6 V a G^VdV a G^] (A.ll) 

With suitable (anti)symmetrization and use of the Bianchi identity \7 a G a a = V a i2, one has, 
for example, 

VaVaGflj = \v {& V {a G m + lc M V ( jG m + l - Cpa C Pa V 2 R (A.12) 

At 6 = the middle term is again expressible in terms of the selfdual part of V a b- The first 
and last terms can be shown to be expressible in terms of the (torsionful) Ricci tensor and 
scalar, as follows: we have the identity expressing the component curvature tensor R a bcd = 
C^gRab-yS + C-ysR abA/ fi in terms of superspace objects [27, 

~ ^C^[C Pl V^G SYt + C ps V { jG 7)7 ] + \c^C p{l C s)a {V 2 R + 2RR) (A.13) 
where the right hand side is evaluated at 9 = 0. In particular 

V(dV (7 G^| e=0 = 2^. ))75 (A.14) 

and, similarly 

V 2 R + 2RR\ e=0 = -\r^ & ^ 5 (A.15) 

Our final ingredient expresses the component torsion and the torsionful connection in terms of 
the axial vector auxiliary field and the torsionless connection J27| , p§| ]: 



Uabc = ^afec(e) - -EabcdV (A. 16) 
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This allows rewriting the component torsionful curvatures in terms of the ordinary curvature 
tensor and additional V a b terms. 

Assembling all the ingredients, we finally obtain the following contributions to the trace and 
/2-current anomalies: 

6 = ^ [{Wa ^ )2 + ( ^ )2] + 3^ [( ^ )2 + ( ^ )2] 

d a Ra = ^[(W a ^) 2 - (W ap . s f\ - ^^[(V aP f - {V ap ?\ (A.17) 

in terms of selfdual and anti-selfdual parts of the curvature tensors and the axial vector field 
strength. This can be rewritten in the form given in ( |4.3| ) and fl4.7|) . 
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